Josephson (001) tilt grain boundary junctions of high temperature superconductors 
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We calculate the Josephson critical current Ic across in-plane (001) tilt grain boundary junctions of 
high temperature superconductors. We solve for the electronic states corresponding to the electron- 
doped cuprates, two slightly different hole-doped cuprates, and an extremely underdoped hole- 
doped cuprate in each half-space, and weakly connect the two half-spaces by either specular or 
random Josephson tunneling. We treat symmetric, straight, and fully asymmetric junctions with s-, 
extended-s, or d^2_y2-wa.ve order parameters. For symmetric junctions with random grain boundary 
tunneling, our results are generally in agreement with the Sigrist-Rice form for ideal junctions that 
has been used to interpret "phase-sensitive" experiments consisting of such in-plane grain boundary 
junctions. For specular grain boundary tunneling across symmetric junctions, our results depend 
upon the Fermi surface topology, but are usually rather consistent with the random facet model of 
Tsuei et al. [Phys. Rev. Lett. 73, 593 (1994)]. Our results for asymmetric junctions of electron- 
doped cuprates are in agreement with the Sigrist-Rice form. However, our results for asymmetric 
junctions of hole-doped cuprates show that the details of the Fermi surface topology and of the 
tunneling processes are both very important, so that the "phase-sensitive" experiments based upon 
in-plane Josephson junctions are less definitive than has generally been thought. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.78.Bz 
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I. INTRODUCTION 

Because of its relation to the mechanism for super- 
conductivity in the high temperature superconducting 
compounds (HTSC), there has long been a raging de- 
bate regarding their orbital symmetry of the super- 
conducting order parameter (OP)i2i2i4 Only "phase- 
sensitive" experiments involving Josephson tunneling 
can distinguish the OP from the non-superconducting 
pseudogap^*^ Although the first in-plane phase- 
sensitive experiment on YBa2Cu307_5 (YBCO) sug- 
gested a dominant s-wave OP^ experiments using 
tricrystal films of YBCO, BizSraCaCuaOg-Ha (Bi2212), 
and Ndi.85Ceo.i5Cu04+5 (NCCO), and tetracrystal films 
of YBCO and La2-^Ce^Cu04- 
that these materials had a dominant 

Qp 1 .7.8.9.10.11.12 

Very different results were obtained for c-axis junc- 
tions. Low temperature T c-axis Josephson junctions 
between Pb and YBCO, Bi2212, and NCCO suggested 
varying amounts of an s-wave OP component i^iiiiiSiiSiii 
More recently, three Bi2212 c-axis twist Josephson junc- 
tion experiments showed that the OP has at least a sub- 
stantial, and possibly a dominant s-wave component for 
T up to the transition temperature y^4*iSii2422i2i Here 
we search for a reason for these qualitative differences^ 

To date, most theoretical treatments of in-plane (001) 
tilt grain boundary (GB) junctions either used the 
Ginzburg-Landau (GL) approachfSS or assumed a circu- 
lar in-plane Fermi surface (PS) cross-section, and that the 
FS-restricted 42_j,2-wave OP cx cos(20fc).^3,24^25,26,27 
Most of those treatments focussed upon the single quasi- 
particle density of states on the surface, rather than the 
critical current across the misalignment grain boundary 



junctions. The only previous treatment that included 
a tight- binding FS was that of Shirai et al., for which 
the FS consisted of very small pockets centered at the 
corners of the first Brillouin zone (BZ)i^ Their chosen 
FS was extremely different from any that have been ex- 
tracted from angle-resolved photoemission spectroscopy 
(ARPES) experiments. Shirai et al. studied both the 
single quasiparticle and supercurrent properties of GB 
junctions constructed from four specific surfacesi^ Here 
we explicitly take account of the tight-binding hole and 
electron-doped FS's observed using ARPES and of the 
surface boundary conditions (BC's) at the interfaces. 

We find that a d^-2_y2-w&Ye OP can be consistent with 
the tricrystal experiments for electron-doped cuprates, 
for which the two-dimensional (2D) FS cross-section is 
nearly isotropic. For hole-doped cuprates, the situation 
is found to be inconclusive, as one FS cross-section is 
consistent with the tricrystal experiments for random 
GB tunneling, but not for specular GB tunneling. An- 
other, slightly different hole-doped FS, on the other hand, 
is consistent with the tricrystal experiments for specu- 
lar tunneling, but not for random tunneling. The FS 
of Shirai et al. is inconsistent with most experiments 
on the hole-doped cuprates, including the tricrystal ex- 
periments, for both specular and random GB tunneling. 
These combined results suggest that defects and small 
changes in the FS topology play essential roles in inter- 
preting the in-plane GB junction experimental resultsi 
An important modification to the tetracrystal experi- 
ment is also warranted.^ 

We let 9i for J = L,Rhe the angles the Xi axes on 
the left (L) and right (R) sides make with the normal to 
a straight GB interface. Fig. 1. GL treatments of the 
Josephson current / for a c?2.2_j^2-wave OP yielded for 
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FIG. 1: (Color online) Plots of /co(^)//co(0) for a ci^2_y2-wave 
OP with straight (dotted green), symmetric (dotted green) 
and asymmetric (solid blue) GB junctions in the ideal GL 
model, and straight (dot-dashed black), symmetric (dashed 
red) and asymmetric (solid blue) GB junctions in the faceted 
GL model. Inset: Sketch of a (001) tilt GB junction. 



ideal (non-reconstructed) and faceted GB's 



7.22 



II = /^rfcos(2^i)cos(2^,j)sin(50), 



/^rfCos[2(^^L+0fl)]sin(,50)/2, 



(1) 
(2) 



where Scj) — (t>L-' 4>r is the OP phase difference across the 
GB. These have the form / = Ico sin((S0), where Ic = \Ico\ 
is the critical current. Ico > and Icq < are the "0- 
junction" and "7r-junction" cases, as I — Ic sin{S(f> + n) 
in the latter. There are three types of junctions we shall 
consider. These are asymmetric, symmetric, and straight 
junctions. We note that other workers have used the 
terminology "mirror" and "parallel" for symmetric and 
straight junctions, respectively.^^ For asymmetric junc- 
tions, Ol ^ e,eR = 0, I'/iO)/r/iO) = cos(26l). For 
symmetric junctions with Ol — Or — 9 and straight junc- 
tions with '6Ih = -9l = 0, /^(6')//d(0) = cos2(26'). For 
faceted symmetric and straight junctions, I^{9)/l[{0) = 
cos(40), 1, respectively. These ideal and faceted GL re- 
sults shown in Fig. 1 are in qualitative agreement with 
our results for electron-doped cuprates with random and 
specular tunneling, respectively, but differ qualitatively 
from our microscopic results for hole-doped cuprates. 



II. PROCEDURE 

To calculate / across an in-plane GB, we assume the 
Fourier transform of the quasiparticle Green function ma- 
trix in the bulk of the ith superconductor is2£ 



iuj + A(fc')Tl + ^(fc»)T3 
'w2 + |A(fc^)|2+e2(feO' 



(3) 



FIG. 2: (Color online) Plots of the four 2D Fermi surfaces 
studied: The nearly isotropic FS3 (thick solid black), the 
holed-doped FS2 (dotted blue)^ the best fit to the ARPES 
FS of Bi2212 (dashed redr^ and the extreme hole-doped FS 
studied by Shirai (thin solid green) i— 



^o(fc*) = -J||[cos(A:>) H- cos(fc^a) 
— t/cos(/c*a) cos(fc*a) — /i]. 



(4) 



where w is a Matsubara frequency, is the wave vec- 
tor on the ith GB side, A(fc') is the respective two- 
dimensional (2D) gap function, a is the tetragonal in- 
plane lattice constant, /i, J||, and v are parameters defin- 
ing the chemical potential, and in-plane bandwidth and 
dispersion, respectively, and the Tj are the Pauli ma- 
trices, g is the usual rank 2 matrix with elements 
G,F,—G^, and that describes both the single quasi- 
particle and mean-field pair excitations. We take ^(fe*) = 
£,o{k^) + J±[l — cos{kls)], where s is the c-axis repeat dis- 
tance. Here we are mainly interested in the 2D limit 
J± = 0. For one tight-binding 2D FS appropriate for 
Bi2212, FS2, we take J|| = 500 meV, = 0.6, and 
1/ = 1.3, and set ^o{k},) = 0.^ A slightly different FS, 
denoted ARPES, that most closely resembles the FS of 
Bi2212 as observed in ARPES experiments, is obtained 
with J|| = 306 meV, /i = 0.675, and = 0.90^28: A nearly 
circular FS, FS3, is obtained with /i = 1.0 and ly = Qm^ 
Finally, we also studied the extreme hole-doped near- 
neighbor tight-binding FS used by Shirai et al^ Taking 
the maximum value of the gap to be 20 meV, this FS is 
governed by the parameters J|| = 201 meV, 1^ = 0, and 
fi = —1.94. These four FS's are pictured in Fig. 2. 

Next, we construct the half-space Green function ma- 
trices on the ith side of a straight (001) tilt GB. We 
let fi and e be unit vectors normal and parallel to the 
GB satisfying e x h = z, as sketched in Fig. 1. Let- 



we set k|| = (fc||, fc*), 



II' 
cos( 



ting /cj*! — ki ■ e and k^^ = ki ■ n, 

— k^± cos 9i — fc^i sin6'i, and ky — n,j^oi±i.yi t n.^^^ 
The lack of translational invariance requires a discrete 
indexing of the lattice planes along n, as in our previous 
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technique for c-axis tunneling:^''. Thus, for integers n, m, 

'"'/"^-^ dk\ 



zA;^ a\ (m — n) 



where fcf 



— (fc^|,/c*) and the a^j_/a — min(cos6'j, sinfl^) 

for 90° > 61, > 0°, = a for 6*, = 0°,90° arc the 
respective lattice plane separations along n. We choose 
the GB interface to be between lattice planes 1 and 
and set n,m > 1, n,m < for i — R,L, respectively. 
Suppressing the and co dependencies, we construct 

each half-space g„m from a combination of the Smn that 
obeys the free surface boundary conditions (BC's) gmo = 
gon = for n, TO > 1 and gmi — gin = for n, to < 0, 
obtaining'^^'^- 



1- 



Q O 

o o 

(110) 



(100) 



1 

x/a 



— G\m~7i\ 
= G\m-n\ 



Gra+n i 
G2—m—n ; 



for m,n > 1, (6) 
for TO, n < 0. (7) 



We then paste together the two half-space layered su- 
perconductors. The supercurrent / across the GB be- 
tween the two half-spaces is given for a general junction 
hopping matrix element Jj{kP^ , kl!) by 



ieT 



E 



Tr 



(to -I- T3) 



x( Jj(fc^,4)Gio(fe^,fc^,A;„c^) 



-J.7(4,fc^)Goi(fc^,fe^,fc.,c^) 



(8) 



where GnmiJ^y^-, kz, is the full Green's function ma- 
trix for the combined two half spaces coupled via tunnel- 
ing across the (001) tilt GB.^° An example of the re- 
sulting interface is pictured in Fig. 3. In principle, it is 
straightforward to evaluate / to all orders in Jj, as was 
done for the case of coherent tunneling across a c-axis 
twist junction/^ As in that case, however, the most rele- 
vant case, especially for GB misalignment angles 9r + 9l 
sufficiently large, is for weak tunneling across the GBn^ 
We set A(k*) = exp(i0,)ReA(k') and define gf,l^ to be 
the (j)q)th matrix element of the rank 2 matrix g„„i with 
(j^i = 0. To leading order in the general tunneling matrix 
element Jj(k/j', k/^), we obtain 



/ = 



11 



Y: |J.(k^,kff)pg}?(kff,..) 



X5^i(kf|,t^) sm{S(t>), 



(9) 



where A is the junction area, e is the electron charge, 
5(t) = (j)L - (l)R, -tt/s < kl < tt/s, -7r/a|| < k'^^ < 7r/a||, 
and the a|| = a^/a^ are the lattice constants along e. We 

note that 5}^ and ggj are the F and pair functions for 
real OP's on opposite sides of the junction, respectively. 



FIG. 3; (Color online) Sketch of a small part of a (llO)i(lOO) 
junction. The bottom rows of the left (110) lattice (open 
red circles) and right (100) (filled blue circles) lattice are 
aligned. An incident quasiparticle can hop along the path 
i, and be specularly reflected along the path r, but tunnels 
non-specularly into the (100) lattice via the two primary ini- 
tial paths ti and t2, as indicated. 



We note that although Eq. ^ is proportional to sin (50 
since it is obtained in the limit of weak GB tunneling, 
our exact expression for the supercurrent, Eq. (|SJ, also 
contains the full temperature dependence of all of the 
harmonics proportional to sin{nd(j)) for n > 1. 

Experimentally for large symmetric tilt angles 6 > 
Oc ~ 5-7°, Ic{e) < /c(0))^-?^i?? justifying our weak 
tunneling assumption, Eq. (Q, although the approx- 
imation is less accurate for 9 < dc- Others working 
on (001) tilt grain boundary junctions of HTSC have 
also made this assumption^SS*^ although they used a 
circular FS and did not necessarily take account of the 
surface boundary conditionsi^ Since the effects induced 
by the GB to the real-space dr^2_y2-wave pairing inter- 
action are formidable, here we only treat the surface 
BC effects, and discuss the effects of self-consistency 
elsewhere. ■^^ We note that a self-consistent calculation 
by Tanaka and Kashiwaya for an OP of mixed symmetry 
suggested no qualitative changes upon the imposition of 
self-consistencyn^ We include the limiting cases of specu- 
lar (or coherent) tunneling, where fc|^ — fe|^, and random 
(or incoherent) tunneling, where fc|^ and fc|^ are indepen- 



dent, writingiii 



R\\2 



\jj{ki:\,^)\ 



(10) 



We remark that the specular (coherent) and random 
(incoherent) tunneling processes are microscopic pro- 
cesses that take place on an atomic scale. For specular 
tunneling to occur, the component of the wave vector of 
the quasiparticle parallel to the junction does not change 
at each atomic site upon tunneling. With random tun- 
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neling, the components of the quasiparticle wave vectors 
parallel to the junction on the two sides of the junction 
are completely independent variables. For example, let 
us consider a tight-binding model of near-neighbor hop- 
ping in a tetragonal lattice, as in Eq. (4) with ly — 0. 
In this model, a quasiparticle far from the junction hops 
from one site to an adjacent site. In a defect-free straight 
(or symmetric) (100) | (100) junction, the lattices on both 
sides of the junction are identical. Thus, a quasiparticle 
can propagate from any site on one side of the junction 
to the nearest neighbor site on the other side without 
changing its direction. A quasiparticle can also backscat- 
ter coherently at the interface, preserving the compo- 
nent of the wave vector parallel to the junction, as in 
other treatmentsi^^^ For instance, when the interface 
is a (110) surface, the surface is theoretically atomically 
flat, allowing a substantial amount of coherent reflection, 
as sketched in Fig. 3. Although with the usual vapor 
phase epitaxy, such (100)|(110) junctions are far from 
perfect ir when they were constructed using the far supe- 
rior liquid phase epitaxy, an excellent fit to the Fraun- 
hofer pattern was obtained by application of a magnetic 
field parallel to the junction, without any evidence of a 
zero-bias conductance peak at the interface^^ at least in 
some cases. Such behavior is consistent with the incoher- 
ent tunneling model of Ambegaokar and Baratoff (AB) 
for an s-wave superconductor. In earlier vapor deposition 
of the somewhat different (110)|(001) GB's, sometimes a 
rough approximation to the Fraunhofcr pattern was ob- 
served, with a peak at low fields, and sometimes a very 
different pattern with a dip at zero field, possibly con- 
sistent with a d-wave OP,^'* was observed4& In both of 
these experiments, high resolution transmission electron 
microscopy (HRTEM) demonstrated the quality of the 
j unctions i^Siffii 

However, at interfaces such as those containing a (130) 
surface. Fig. 4, the amount of coherent reflection off 
the surface is likely to depend strongly upon the inci- 
dent quasiparticle wave vector direction, as the surface 
is a periodic, but uneven zigzag. To date, such junctions 
have never been prepared with the atomic precision of 
the above two casesi^ 

On the other hand, for an asymmetric (100)|(110) junc- 
tion. Fig. 3, the ratio of the parallel lattice spacings on 
the two sides of the junction is v^, an irrational number. 
Although specular reflection at the interface by a quasi- 
particle is entirely possible, in this case, even if the junc- 
tion were atomically perfect, the proportion of atomic 
sites at which coherent tunneling across the junction can 
occur is vanishingiy small. That is, in tunneling to the 
nearest site across the junction, nearly all quasiparticles 
will change the component of their direction parallel to 
the junction. One example of this is sketched in Fig. 3. 
When one includes the dominant junction tunneling from 
all sites along the junction, the tunneling component of 
the quasiparticle wave vector parallel to the junction is 
effectively random. 

In general, one expects that atomically perfect straight 
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FIG. 4: (Color online) Sketch of the atomic sites (filled red 
circles) of a tetragonal plane near a (130) surface (the x axis), 
where a' — vTOa. The solid black and dashed blue arrows in- 
dicate selected incident and reflected quasiparticle paths in- 
volving only near-neighbor and next-nearest neighbor hop- 
ping, respectively. 



junctions can exhibit specular, or coherent tunneling. 
Defects would lead to a component of the tunneling that 
is random, or incoherent. Asymmetric junctions and 
symmetric junctions that cannot also be described as 
straight junctions (such as the (100)|(100) junction) in- 
trinsically lead to predominantly random, or incoherent 
tunneling. Thus, atomically perfect (jjqO)\{100) asym- 
metric or {pqO)\(pqO) symmetric junctions for p + q ^ 1 
(with 9l — Or) are expected to exhibit primarily random, 
or incoherent microscopic tunneling processes. Facets, 
atomic disorder, stoichiometry variations, etc., at the 
junction contribute further to the randomness of the tun- 
nehng processesM-^ii^24iM4^ 

This microscopic (atomic scale) coherence or incoher- 
ence of the tunneling process is completely unrelated to 
the macroscopic coherence of the superconducting wave 
function, which can occur over the entire area of the junc- 
tion, provided that it is prepared sufficiently uniformly. 
For example, the AB microscopic model of Josephson 
tunneling, '^^ which averages independently over the par- 
allel components of the quasiparticle wave vector on each 
side of the junction, assumes uniformly random, or inco- 
herent, microscopic quasiparticle tunneling only. How- 
ever, that model leads to macroscopic coherence of the 
superconducting wave function on each side of the junc- 
tion. Note that although most experiments are consis- 
tent with Bi2212 having strongly incoherent intrinsic c- 
axis tunneling between the double layers, the product of 
the critical current and the normal state resistance IcRn 
is about 20% of that expected for a pure s-wave super- 
conductor with perfect interfaces, the AB value4i3^ This 
occurs both for intrinsic and for c-axis twist junctions 
several fj,m in cross-section4iiSii2i2i 
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FIG. 5: (Color online) Plots of logj^,, \Ic{6)/Io\ for asymmetric 
junctions with Or = 0, 8l ~ 0, for s- (thin blue), extended-s- 
(intermediate thickness red), and d^2_y2-wave (thick black) 
superconductors, for specular (solid) and random (dashed) 
GB tunneling and the tight-binding hole-doped FS2. The 
d-wave curves satisfy 7co(90° — 6) = —IcoiO), but the s- and 
extended- s- wave curves satisfy Jco(90° — 0) = Ico{6). See text. 



III. RESULTS 

In Fig. 5, we present our results obtained 
from Eq. © for the low-T across asymmet- 
ric junctions, 9ji = 0, 0^ = 9, including the 
tight-binding FS2 and the surface BC effects. In 
all of our calculations, we assume the d^2_y2- and 
extended-s-wave OP's have the bulk real-space pairing 
forms Ao(T)[cos(fc*a) - cos(fc^a)] and Ao(T)| cos(fc*a) - 
cos(fc^a)|, respectively. The ordinary- s- wave OP is 
Ao(r), independent of k'M We take Ao(0) = 10 meV, 
which implies a maximum d-wave OP of 20 meV. In 
Fig. 4, we plotted log^o |/c(^')/-?^o|, where 4 = |/co| 
and lo = 4e|J7|V[a3s(0.1J||)2],4e|J}|V[a4s2(o.iJ||)2], 

respectively!^ The solid (dashed) curves are the results 
for specular (random) tunneling, respectively. The spec- 
ular d-wave GB junction behaves as a 7r-j unction only for 
14.9° < 6* < 45° and 75.1° < 9 < 90°, qualitatively dif- 
ferent from the GL results pictured in Fig. l.^^ Note that 
the d-wave Ic for random GB tunneling is small but non- 
vanishing due to the surface BC's. For this and all other 
cases of random tunneling across GB's between (1^2 _y2- 
wave superconductors, neglect of the surface BC's would 
lead necessarily to Ic = 0. In the 2D limit with random 
GB tunneling, however, the asymmetric GB results are 
qualitatively in agreement with the GL results, except 
for the overall magnitude of Ic, which is greatly reduced 
from that obtained for specular GB tunneling. 

We also investigated the role of J± for the d-wave case. 
For random GB tunneling, as J± increases to 40 meV, the 
broad peak in Ico{9) becomes much flatter, and for much 



FIG. 6: (Color online) Plots of logj^Q \Ic{0)/Io\ for symmetric 
junctions with 9r — 6l = 9, for s- (thin blue), extended- s- 
(intermediate thickness red), and d^2_y2-wave (thick black) 
OP's for specular (solid) and random (dashed) GB tunneling 
with the tight-binding, hole-doped FS2. Results for straight 
junctions with Or = —9l = d and a d^2_y2-wave OP for 
specular (dotted green) and random (also dashed black) GB 
tunneling are also shown. All curves satisfy /co(90° — 9) = 
Ico{9). 



larger Jj^, Icf){9) non-monotonically approaches the GL 
result oc cos(2^?). For specular tunneling, = 10 meV 
yields d-wave results nearly identical to those pictured for 
Jj^ = 0. Increasing J±^ beyond J_lo ~40 meV eliminates 
the regions within Q < 9 < 45° of 7r-junctions, so that 
the overall region of 7r-junctions is for 45° < 9 < 90°, 
as for the GL case. Hence, for J± < J±o, the specular, 
tight-binding asymmetric GB d-wave results with FS2 
are qualitatively different than the GL ones. 

In Fig. 6, we present the analogous results for sym- 
metric junctions with FS2. We also plotted our re- 
sults for straight GB junctions with a dx2-y2-wa,ve OP. 
For random GB tunneling, the straight GB junction re- 
sults are identical to the dashed black curve for random 
GB tunneling with symmetric junctions. For specular 
GB tunneling, the straight junction results are shown 
as the dotted green curve. Note that for specular tun- 
neling, straight and symmetric junctions have opposite 
Ic values at 6* = 45°. For symmetric GB junctions, 
the specular d-wave junction behaves as a 7r-junction for 
39.4° < 9 < 50.6°, but behaves as a 0-junction other- 
wise. For random tunneling, the symmetric d-wave GB 
junction behaves similarly to the ideal GL model, Eq. 

as both have /c(45°) = but never behave as tt- 
junctions. However, the d-wave GL facet model, Eq. (2), 
has /J (45°) — — /J(0°), which is qualitatively different^ 
We note from Figs. 5, 6 that for small angle {0 < 9 < 5°) 
specular GB junctions, the extended-s and dj.2_j,2-wave 
OP's lead to indistinguishable Ic{9)- However, for ran- 
dom GB tunneling, the d-wave /c(^) is much smaller than 
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FIG. 7: (Color online) Plots of logj^,, jJc(6)//o| for symmetric 
junctions with specular (thick solid red) and random (thick 
dashed black) tunneling, and asymmetric junctions with spec- 
ular (thin solid blue) and random (thin dashed purple) tunnel- 
ing, for a d^2_j^2-wave OP with the ARPES FS.^^ The straight 
and symmetric GB junctions satisfy 7c (90° ~ 9) — lc{0), but 
the asymmetric GB junctions satisfy 7c(90° — 6) = —lc{0). 



the extended-s-wave Id^), and would vanish in the ab- 
sence of the surface BC's. Except for the broad peak and 
the greatly reduced magnitude of Ic, the symmetric GB 
junction across d^2_y2--wave superconductors with FS2 
and random GB tunneling is qualitatively in agreement 
with the GL result of Fig. Im^ For specular GB tunnel- 
ing, however, our microscopic dx2-y2 results differ quali- 
tatively from the GL formsiSS 

We remark that the ordinary-s- and extended-s-wave 
OP's lead to /c(0) values that are comparable to each 
other for both GB tunneling mechanisms. However, for 
the da.2_y2-wave OP, random GB tunneling with FS2 
leads to a suppression of /c(0) by a factor of 461 from 
that of the ordinary-s- wave OP. This suggests that /c(0) 
for the dj,2_j,2-wave OP with weak, random GB tunneling 
is unlikely to be consistent with the large values observed 
experimentally across low-angle (001) tilt junctions.^ The 
broad peak in /co(0) for the d-wave OP with random GB 
tunneling for both asymmetric and symmetric junctions 
is also inconsistent with most experimental data. On the 
other hand, the d-wave OP with specular GB tunneling 
across asymmetric junctions does show a rapid decrease 
in lco{0) / IcoiO) with increasing 6, in agreement with low- 
angle experiments,^ but vanishes at 14.9° and 39.4° for 
FS2, respectively, unlike the experiments.^ 

In Figs. 7-9, we compare our log^Q |/c(^)/-fo| results 
for dx2-y2-'wave superconductors across symmetric and 
asymmetric GB junctions with three different FS's. For 
each FS, Ic results for straight and symmetric GB junc- 
tions with random GB tunneling are identical, and for 
specular GB tunneling, the /c(45°)'s are opposite in sign 
and equal in magnitude. In Figs. 7-9, the results for sym- 



metric GB junctions with specular (random) GB tunnel- 
ing are shown as the thick sohd red (dashed black) curves, 
respectively, and results for asymmetric GB junctions 
with specular (random) GB tunneling are shown as the 
thin solid blue (dashed purple) curves, respectively. Re- 
sults for straight and symmetric GB junctions with ran- 
dom tunneling are identical, and the results for straight 
GB junctions with specular GB tunneling are shown as 
the dotted green curves. 

In Fig. 7 we plotted log^Q \Ic{0)/Io\ versus 9 for sym- 
metric and asymmetric GB junctions for a dx2-y2-wave 
superconductor with the ARPES FS. Note that for sym- 
metric junctions, a 7r-junction occurs with specular GB 
tunneling for 40.3° < 9 < 49.7°, very similar to the re- 
gion obtained for FS2 pictured in Fig. 5. In addition, 
for random GB tunneling across symmetric junctions, no 
TT-junctions are obtained. However, for asymmetric junc- 
tions, the results differ substantially from those obtained 
for a dx2^y2-'W8ive superconductor with FS2, pictured in 
Fig. 4. In that (FS2) case, the asymmetric GB junc- 
tion was a TT-junction for specular GB tunneling with 
14.9° < 9 < 4:5° and 75.1° < 6* < 90°, but was a 0- 
junction otherwise. For random GB tunneling, it was 
always a 0-junction for < < 45° and a 7r-junction for 
45° < 6* < 90°. For this slightly different (see Fig. 2) 
ARPES FS, however, nearly the opposite behavior was 
found. Asymmetric GB junctions with the ARPES FS 
with specular GB tunneling are 0-junctions for < 6* < 
45° and 7r-junctions for 45° < 9 < 90°, but for random 
tunneling, a 7r-junction appears for 30.6° < 9 < 45° (and 
also for 59.4° < 9 < 90°), and a 0-junction otherwise. 
Hence, these slightly different hole-doped FS's lead to 
qualitatively different Ic behavior for d^2_y2--wave super- 
conductors. For s- and extended-s-wave superconduc- 
tors, on the other hand, such strong FS dependence of 
Ic{9) does not occur. Although for clarity in present- 
ing the d-wave results, we have not shown our s- and 
extended-s-wave data for the ARPES FS, our calcula- 
tions confirm the close similarity of those results to the 
analogous ones pictured in Figs. 5 and 6, confirming the 
reliability of our results. 

With the ARPES FS, asymmetric GB junctions with 
specular tunneling and for symmetric GB junctions with 
random tunneling are qualitatively consistent with the 
ideal GL model, and symmetric GB junctions with 
specular tunneling are qualitatively consistent with the 
GL facet model, ^ except that the position of Ic{9) — 
has shifted from 22.5° to 40.3°. However, asymmetric 
GB junctions with random tunneling and the ARPES 
FS are completely unlike either GL model. 

Next, we calculated Ico{9) at low T across GB's with 
the nearly circular FS3, appropriate for the electron- 
doped cuprates such as NCCO, and presented our re- 
sults in Fig. 8. For asymmetric GB's, our d-wave re- 
sults agree with the GL model, /co(^)/^co(0) ~ cos(20), 
for both specular and random GB tunneling, but for the 
latter, /co(0) is greatly reduced in magnitude for a d- 
wave OP. Our results for symmetric or straight GB's 
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FIG. 8: (Color online) Plots of logj^Q \Ic{0)/Io\ for symmetric 
junctions with specular (thick solid red) and random (thick 
dashed black) tunneling, and asymmetric junctions with spec- 
ular (thin solid blue) and random (thin dashed purple) tunnel- 
ing, for a d^2_y2-wave OP with the nearly isotropic, "electron- 
doped" FS3. Also shown is the specular d-wave curve (dot- 
ted green) for straight junctions. The straight and symmetric 
junctions satisfy Ic{90° ~9) = lc{0), but the asymmetric junc- 
tions satisfy //(90° -6) = -Ic{e). 



with random GB tunneling coincide with each other, 
and Ico{9) / lco{0) ~ cos^(26'), as in the GL model for 
ideal symmetric GB junctions (Fig. 1). For specular 
GB tunneling across symmetric GB junctions, the re- 
sults are qualitatively similar to the GL facet model, 
Eq. (2) except that the 9 value for the crossover from 
0- to TT-junction behavior has shifted from 22.5° to 29.8°. 
With that modification, all of our d-wave results for the 
electron-doped FS3 are qualitatively in agreement with 
the GL facet model for specular GB tunneling^ and with 
the GL ideal model for random GB tunneling between 
dx2_y2-wave superconductors.^^ This conclusion is some- 
what unexpected, as one naively might have expected 
the specular and random results with FS3 to agree with 
the ideal and facet GL models, respectively. Our low-T 
symmetric junction results are also qualitatively in agree- 
ment with those obtained at low T for weak tunneling 
with a perfectly circular FS and a FS-restricted d-wave 
OP with different GB roughness, and our straight GB 
junction results with random GB tunneling are quali- 
tatively consistent with previous straight, dirty-A^, SNS 
junction resultsr^ 

In order to compare the results obtained using our 
technique of imposing the surface BC's upon the Green 
function matrix to that obtained using different tech- 
niques by others, we performed explicit calculations for 
these purposes. According to the formulas in their pa- 
per, Shirai et al. used the FS pictured by the small green 
pockets near the corners of the first BZ pictured in Fig. 

Those authors constructed half-spaces exhibiting 
(100), (130), (120), and (110) surfaces, and pasted them 



together with specular tunneling in either the straight or 
symmetric configurations (which they denoted as paral- 
lel and mirror, respectively). Their discrete quantization 
axis was always along the x direction parallel to the CuO 
bonds, and hence only normal to the surface for the (100) 
case. The number of rows in the y direction varied from 
1 for the (100) and (110) surfaces to 3 for the (130) sur- 
face. Unlike our quantization axis that was defined to 
be parallel and normal to the surface as pictured in Figs. 
3 and 4, upon Fourier transformation in the y direction, 
their non-normal quantization lattice basis set could lead 
to unintended hopping across the GB with the strong, 
intrinsic tight-binding near-neighbor hopping matrix el- 
ement. 

In addition, Shirai et al. included tunneling from a 
few rows other than the surface row of idealized atomic 
sites, by assuming that the specular GB tunneling ma- 
trix elements decreased inversely instead of exponen- 
tially with the tunneling distance, taking the maxi- 
mum hopping strength to be Q.05t, or 6.2% of the ef- 
fective near-neighbor hopping strength, 0.803i.^^ Never- 
theless, they calculated the single quasiparticle densities 
of states and Ic across the straight and symmetric GB 
junctions with the above form of specular tunneling, for 
9 = 0°, 18.4°, 26.5°, and 45°, respectively. Some of their 
low-T limits of Ic can be compared with ours using their 
unusual FSi^S Although they did not mention it explic- 
itly, Shirai et al. found the straight 9 = 0° and 26.5°, 
junctions to be 0-junctions. However, they found strong 
zero-bias conductance peaks at all junction sites for the 
straight 45° junction, and at two of the three distinct 
junction sites for their 18.4° junction, suggesting that 
both of these junctions behave as 7r-junctions. For the 
symmetric junctions, they found the low-T limits of the 
symmetric 18.4° and 45° junctions to be 7r-junctions, but 
the low-T limits of the 0° and 26.5° symmetric GB junc- 
tions to be 0-junctions.2£ Of these, Ic(T) for the 18.4° 
symmetric junction changed sign at a finite T value. 

Our low-T results for log^Q \Ic{9)/Io\ versus 9 for sym- 
metric, asymmetric, and straight junctions with the Shi- 
rai FS are pictured in Fig. 9. In all cases, even for 
9 = 0, the extremely small overlap of the Fermi surfaces 
for specular GB tunneling leads to an extremely small 
Ic- Our results for specular GB tunneling with sym- 
metric junctions are shown as the thick solid red curves 
with three alternating regions of tt- and 0-junctions for 
< 6* < 45°. In the domain 0° < 6* < 45°, 7r-junctions 
are found in the three regions 1.57° < 9 < 7.19°, 
9.85° < 9 < 20.89°, and 31.52° < 9 < 45°. The be- 
havior of the domain 45° < 9 < 90° is obtained from 
/c(90° — 9) = Ic{9)- Random tunneling across either 
straight or symmetric GB junctions leads to just a 0- 
junction for < ^? < 90°, although there are strong 
dips in Ic at 9 fi: 13.6° and 76.4°) (not pictured) and 
/c(45°) = 0, as indicated by the thick dashed black 
curve. For asymmetric GB junctions with specular GB 
tunneling, the thin solid blue curve shows a 7r-junction 
for 17.04° < 61 < 45° (and from 72.96° < 9 < 90°), 



8 




FIG. 9: (Color online) Plots of logj^Q \Ic{0)/Io\ for symmetric 
junctions with specular (thick solid red) and random (thick 
dashed black) tunneling, and asymmetric junctions with spec- 
ular (thin solid blue) and random (thin dashed purple) tun- 
neling, for a d^'2_y2-wsye OP with the extreme hole-doped 
FS of Shirai et alm^ Also shown is the specular d-wave curve 
(dotted green) for straight junctions. The straight and sym- 
metric junctions satisfy 7c(90° — 9)= lc{0), but the asymmet- 
ric junctions satisfy 7c(90° — 8) — —Idfi). The thick dashed 
black curve represents 0-j unctions only. 



FIG. 10: (Color online) Plots of Ico{9)/ho{Q) for the d^2_y2- 
wave OP with specular tunneling, with (solid) and without 
(dashed) the surface boundary conditions, for symmetric GB 
junctions with the nearly circular FS3 (thick black) and the 
tight-binding FS2 (thin red) and , and for asymmetric GB 
junctions with FS2 (intermediate thickness blue). The sym- 
metric curves satisfy Jc (90° — 9) = Ic{9), and the antisym- 
metric curves satisfy Jc (90° — 9) — —Ic{9). 



and the thin dashed purple curve for random GB tun- 
neling indicates a 7r-junction for 13.34° < 9 < 45° (and 
from 76.66 < 9 < 90°). In addition, the dotted green 
curve represents straight junctions with specular tunnel- 
ing, which are always 0-junctions, as in the ideal and 
faceted GL models pictured in Fig. 1, and in all other 
FS's we studied (e. g.. Figs. 6 and 8). 

Some of our results and those of Shirai can be directly 
compared. Both procedures lead to T = 7r-junctions 
for symmetric GB junctions with specular GB tunnel- 
ing with 9 = 18.4°, 45°, but 0-junctions at 6* = 0°,26.5°. 
As for straight junctions, our results for all four FS's we 
studied always lead to 0-junctions, as for both the ideal 
and faceted GL models. For straight junctions at the 
two angles 0° and 26.5°, Shirai et al. also obtained 0- 
junctions for all T < Tc- However, as mentioned above, 
Shirai et al. appear to have obtained 7r-junctions for 
the straight junctions with 9 — 18.4° and 45°. In par- 
ticular, we find their apparent 7r-junctions for the two 
straight junctions, 9 = 18.4°, 45° to be difficult to un- 
derstand. Suppose Shirai et al. were to increase their 
dominant junction coupling strength from 0.05t to their 
bulk hopping strength, G.803t. Then, the junction would 
be invisible to the quasiparticles, and one would certainly 
not expect any surface effects, such as in the straight 0° 
junction they studied?^ Although Shirai et al. attributed 
their peculiar results to the oddness of the number of site 
rows in their minicell, we suspect that their results may 
be artifacts of their peculiar minicell basis, which was not 



defined to be parallel and perpendicular to the GB, as in 
Fig. 4^ 

Finally, in Fig. 10 we compare our specular 
IcoiS) / lco{0) d-wave results with the surface BC's with 
those obtained from the bulk Green functions. We com- 
pare those results obtained with the tight-binding FS2 for 
asymmetric (Fig. 5) and symmetric (Fig. 6) junctions, 
and also our results for symmetric junctions obtained 
with the nearly circular FS3 (Fig. 8). For both sym- 
metric GB cases, there are no regions in 0° < < 45° of 
TT-j unctions for the bulk calculations. Including the sur- 
face effects, these symmetric GB cases have 7r-j unctions 
for 39.4° < 9 < 50.6° and 29.2° < 9 < 55.8° for 
FS2 and FS3, respectively. For asymmetric junctions 
with FS2, the bulk calculation leads to weak 7r-junctions 
for 23.7° < 9 < 28.3° with a negative slope at 45°, 
yielding additional 7r-junctions for 45° < 9 < 61.7° 
and 66.3° < 6* < 90°. The surface BC's greatly en- 
hance the magnitude and range of such 7r-junctions to 
14.9° < 6* < 45° and 75.1° < < 90°, and change the 
sign of the slope at 45°. Although not pictured in Fig. 
10, the surface BC's play an essential role in the case of 
random tunneling across GB's between d-wave supercon- 
ductors. In the absence of the surface BC's, /c = in 
such cases. Thus, we conclude that the surface BC effects 
cannot be ignored, especially with regard to dx2^y2-wsLve 
superconductors. 
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IV. DISCUSSION 

Tricrystal experiments were performed for a variety of 
cuprate superconductorsi^iSiS Tricrystal (a) consists of 
one asymmetric 0°/60° junction and two 0°/30° junc- 
tions. Tricrystal (b) consists of two 0°/71.6° junctions 
and one symmetric 18.4°/18.4° junction. Tricrystal (c) 
consists of two asymmetric 0°/73° junctions and an 
asymmetric 2°/32° junction»i*i2i It was claimed that an 
odd number of 7r-junctions was observed only for tricrys- 
tal (a).^'^" In the case of Bi2212, only results for tricrystal 
(a) were showniSiS 

First of all, our results for the electron-doped FSB with 
either specular or random GB tunneling are in complete 
agreement with a predominant c?3,2_y2-wave OP as deter- 
mined from the tricrystal experiments on NCCO^ That 
is, one would only expect an odd number of 7r-j unctions 
only for tricrystal (a) , for both random and specular GB 
tunneling. 

However, for hole-doped cuprates, the agreement with 
the tricrystal experiments is less robust. Assuming only 
a (i^2_j,2-wave OPtr our FS2 results pictured in Figs. 3, 
4 show that tricrystal (a) would have two 7r-junctions 
if the tunneling were specular, and one 7r-junction if it 
were random. Tricrystal (b) should have no 7r-junctions, 
as claimed. Tricrystal (c), if the 2°/32° junction can be 
approximated by an asymmetric 0°/34° junction, would 
have one 7r-junction for specular tunneling, and two tt- 
junctions for random tunneling. Hence, our FS2 results 
are only consistent with the observation claims if the tun- 
neling is non-specular. 

On the other hand, the ARPES FS, while only slightly 
different from FS2, leads to very different conclusions 
with regard to the agreement with the tricrystal exper- 
iments. For specular GB tunneling, the ARPES FS is 
in complete agreement with the observations of the three 
tricrystal experiments. However, for random GB tunnel- 
ing, tricrystal (a) is expected to have a very small overall 
/c, as /c ~ for all three asymmetric GB junctions in that 
tricrystal. Although tricrystal (b) would be expected to 
exhibit an even number of 7r-junctions, as claimed to be 
observed, tricrystal (c) would be expected to exhibit one 
TT-junction for random GB tunneling with the ARPES 
FS. Hence, the ARPES FS can only be fully consistent 
with the tricrystal experiments on hole-doped cuprates 
for non-random GB tunneling. 

One interpretation is that if the quasiparticle tunnel- 
ing across the relevant grain boundaries is random and 
the OP were pure d-wave, then the FS would be expected 
to be more like FS2. In the unlikely event that the grain 
boundary tunneling were coherent, then a pure d-wave 
OP would have a FS more like the ARPES FS. Such a 
scenario might be used to argue that the FS of Bi2212, 
usually thought to have the ARPES FS form, is actually 
closer to FS2. However, we caution that FS2 and the 
ARPES FS are one-band approximations to the FS of 
Bi2212, which contains two Cu02 layers thought by many 
to be responsible for the superconductivity, and that no 



tricrystal (c) experiments were performed on Bi2212. In 
addition, this scenario ignores the c-axis twist experi- 
ments, which indicate that essentially all of the c-axis 
tunneling is in the s-wave channel,* If those three exper- 
iments are correct, then Bi2212 contains an OP with an 
s-wave component that is at least 20% of the total. Al- 
though not appreciated by Tsuei et al.^^ in the presence 
of nanoscale disorder that has recently been observed in 
STM experiments^^ a mixed d+s OP is indeed allowedi^ 
Thus, it is possible that Bi2212 could have a dominant 
d-wave OP and still have all of the intrinsic c-axis tunnel- 
ing be in the s-wave OP channel. If this scenario proves 
to be true, then further calculations with a mixed s- and 
d-wave OP would be warrantedi2£ 

With regard to the most-studied material, YBCO, in 
addition to the two Cu02 layers, there is also a con- 
ducting CuO chain layer, which is responsible for the 
orthorhombicity and the strong twinning that occurs in 
thin films. The strong twinning makes (100) and (010) 
films indistinguishable, which can complicate the anal- 
ysis. Second, unless strongly underdoped, YBCO ap- 
pears to be more three-dimensional than are most other 
cuprates, so that 7^ might be large enough to cause 
an important modification to our results. Third, it is 
clear from many other experiments that YBCO also has a 
substantial s-wave component ^^ii^ which has often been 
ignored by many workers In each twin domain, a pur- 
ported dominant d-wave OP would have either a d -t- s or 
a d — s OP mix, and if the twin areas were equal, then 
one might expect the s-wave OP component to make 
a negligible contribution to the measured critical cur- 
rents. However, in c-axis YBCO/Pb Joscphson junctions 
made with highly twinned YBCO thin films, quantita- 
tive agreement to the Fraunhofer pattern was obtained, 
although IcRn was considerably less than the AB valuei-^ 

Thus, the different results of FS2 and the ARPES 
FS indicate that the tricrystal experiments are less ro- 
bust determinations of the symmetry of the order pa- 
rameter in hole-doped cuprates than has been gener- 
ally thought. Small variations in the FS topology and 
the non-specularity of the GB tunneling processes play 
a crucial role that cannot be anticipated a priori. To 
the extent that defects play a role in the tunneling pro- 
cesses, it might be possible to understand the experi- 
mental observations with a substantial s-wave OP com- 
ponent. Especially since the (001) tilt GB junctions are 
known to be increasingly underdoped with increasing GB 
misalignment it is possible that 7r-junctions aris- 
ing from the local magnetism of the Cu^ spins formed 
from oxygen loss at the GB junctions could be important. 

All of our calculations using the Shirai FS are in com- 
plete disagreement with both the tricrystal and tetracrys- 
tal experiments, regardless of the degree of randomness in 
the GB tunneling matrix elements. In addition, the Shi- 
rai et al. results for specular symmetric 18.4°/18.4° GB's 
are in direct conflict with all experiments using tricrystal 
(b)i^i^ Taking all of our results for the various FS's into 
account, it appears that there is a topological crossover 
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relevant for the agreement of the critical currents with 
the GL results with full agreement for the electron- 
doped FS's centered about the F point in the first BZ, 
and complete disagreement with the extreme Shirai PS. 
The cross-over appears to occur for FS's intermediate 
to the ARPES PS and FS2. Hence, for the hole-doped 
cuprates, not only are the details of the tunneling impor- 
tant, but the topological structure of the hole-doped PS 
is also very relevant, and minute details can change the 
interpretation. 

Tetracrystals consisting of highly twinned thin films of 
the hole-doped YBCO and electron-doped LCCO con- 
taining a symmetric 0°/0° junction and a 45°/ ± 45° 
junction, plus two symmetric 22.5°/22.5° junctions were 
madeiiiJ^ For comparison, they also made bicrystals 
with two symmetric 22.5°/22.5° junctions. We also cal- 
culated the straight GB junction case 9ji = —9l = 
for a dx2-y2-wa.ve OP, and our results for FS2 and PS3 
are shown by the green dotted lines in Pigs. 6 and 8, 
and at 45°, are opposite in sign to the respective sym- 
metric junction cases with specular GB tunneling. If the 
tunneling were non-specular, as required for the d-wave 
interpretation with PS2 of the tricrystal experiments, 
then the 45°/ ± 45° junction would have the smallest 
(or vanishing, for random tunneling) Ic value for a d- 
wave OP. If instead the tunneling were specular, as re- 
quired for the ARPES PS c?-wave interpretation of the 
tricrystal experiments on hole-doped cuprates, then that 
junction would allow a small but significant Ic value, but 
would behave either as a tt- or 0-junction, respectively, 
for a d-wave OP, depending upon the particular distri- 
bution of twin domains present at the junction. Thus, 
this tctracrystal would have an indeterminate number 
of TT-junctions for specular tunneling and a d-wave OP, 
but no supercurrent for random tunneling. Similar argu- 
ments apply to the LCCO case. Por a highly twinned 
thin film junction, both specular and random tunnel- 
ing across the 45°/ ± 45° GB junction would lead to 
/c ~ for a d^2_y2-wave OP, contrary to the claims in 
the experimentjii The other two junctions, symmetric 
22.5°/22.5° junctions, would not be 7r-junctions, unless 
the tunneling were specular and the PS were different 
from any that we studied. Hence, the interpretation put 
forward by those authors, which assumed the 45°/ ±45° 
junctions were straight, should be reexaminedfii*i^ as the 
TT-junction could arise extrinsically from defects in this 
very imperfect 45°/ ±45° junction. ^■^^■■^^'42, We urge that 
the experiment be redesigned in the form of a tricrystal, 
omitting the 45°/ ± 45° junction. 

Recently, there have been zigzag in-plane YBCO/Nb 
experiments In these experiments, a YBCO film is 
angularly cut from above with a ramping angle ip « 
15 — 20° into a regular zigzag pattern of nominally ex- 
posed (100) and (010) edges, each edge of equal length 
£ 10 — 40/im. After deposition of a thick layer of Au, 
a thicker layer of Nb was applied. This experiment is a 
variation of the earlier 9 = 90° YBCO/Pb superconduct- 
ing quantum interference device (SQUID) experiment of 



Mathai et al.^ which contained YBCO/Pb junctions 
with a thin Ag interstitial layer along the ramped nom- 
inally (100) and (010) edges of YBCO, respectively. In 
those experiments, one of the junctions functioned as a 
0-junction, and the other functioned as a vr-junctioniiS 
The beauty of the recent experiment is that the zigzag 
array of alternating 0-/7r-j unctions was observed by a 
SQUID microscope to exhibit an ordered array of frac- 
tional vortices, the magnitude of which was within the 
experimental error of <&o/2, where 4>o is the flux quan- 
tum. Although no microscopic analysis of the edges of 
the YBCO was presented, the ion milling technique em- 
ployed necessarily leaves a rough surface on a scale of at 
least 1-5 nm, so that the tunneling across each of these 
edges is most likely incoherent. The thick Au layer could 
modify the overall tunneling process, however. The much 
longer length scales Aj of the Josephson vortices, which 
is only 1-2 orders of magnitude smaller than the edge di- 
mensions, allows the vortices to weakly interact, forming 
a regular array. These experiments are suggestive of a 
dominant d-wave OP. 

However, we recall that Gim et al. studied YBCO/Pb 
squid's prepared with a much wider selection of SQUID 
angles 9^ In those experiments, which did not include 
9 — 90°, the results were inconsistent with either a pre- 
dominant d- or s-wave OP component}^ but were in- 
stead consistent with a p-wave polar state. Apparently, 
the preparation of a 0-junction or a 7r-junction depends 
upon the ramping angle Lp of the bombarding ions, so 
that many pairs of SQUID's with angles 9 and 9 + 180° 
contained one 0-junction and one 7r-junction^ Although 
the nice zigzag experiments do demonstrate the repro- 
ducibility of the 0-junction and 7r-junction preparation 
at fixed If and 9 — 90° with YBCO/Nb junctions, those 
authors did not report the results of 9 and ip variation 
studiesiLia 

We are presently engaged in calculating the single 
quasiparticle and pair excitations at low T , both of which 
are contained in the Green function matrix given by Eqs. 
© and Q . Por brevity we have here only presented the 
Josephson critical current across the grain boundaries, 
the quasiparticle densities of states at the sites on the 
rows adjacent to or near to the junction can be evalu- 
ated from the imaginary part of the single particle g}^^ 
for n = 0, 1, . . . Although the most commonly used pro- 
cedure for obtaining the zero-bias conductance peak on 
a surface is to use the Bogoliubov-deGennes technique^^ 
our Green function technique with the appropriate sur- 
face boundary conditions is equally valid, ^^'^^ but is con- 
siderably more compact and powerful. With our Green 
function technique, we are able to evaluate the critical 
current at arbitrary grain boundary angles, and arc not 
restricted to specific interfaces such as {pqQ)\{p' q' Q) for 
integral p,q,p',q'. It is also much easier to generalize 
the tunneling forms, although we have limited our treat- 
ment to the most commonly studied specular and random 
cases. Although it would be rather easy to generalize the 
tunneling matrix elements to forms intermediate to these 
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limiting cases for symmetric and straight junctionS)2i for 
asymmetric junctions, such a generalization is more difS- 
cult. We have compared our results with those obtained 
by Shirai et al. using a purely numerical technique, and 
some of their numerical results at low T are in qualitative 
agreement with our results obtained from Eq. 0. Our 
results with the nearly circular FS3 are in fairly good 
agreement with those obtained with the GL model and 
with microscopic calculations using a circular FSi22i2^ 
The main drawback with the Green function technique is 
the procedure for making the calculations self-consistent. 
However, as noted by Tanaka and Kashiwaya,"*^ mak- 
ing the results self-consistent does not usually result in 
qualitative changes. We are presently engaged in self- 
consistent calculations for various grain boundaries of the 
{pqO)\{p'q'0) types, as was done by Shirai et al.'^^ but we 
are using the near-neighbor and next-nearest neighbor 
single quasiparticle matrix elements appropriate for ei- 
ther the ARPES or FS2 FS's on Bi2212, and wiU present 
our results elsewherei^ We will also present the quasipar- 
ticle densities of states for these more realistic FS models 
elsewhere 

Finally, we remark that recent YBCO thin film struc- 
tures have been studied using (001)|(103) thin film 
boundaries with 45° in-plane misalignment angles, and 
the results interpreted using the Sigrist-Ricc formulai^2i^ 
We have not yet analyzed such grain boundaries, be- 
cause the Fermi surface of YBCO is particularly com- 
plicated, and our FS2 and ARPES FS's are likely to 
be inadequate to do so reliably. Nevertheless, from 
the results we obtained with these two FS's, it is ev- 
ident that the results obtained from such devices, the 
grain boundaries of which have not been analyzed with 
atomic scale precision as was done with the Bi2212 c-axis 

twist experiments^ might also have a more conventional 
explanation.5ii52i53,54,55,56 

It has been suggested that the experimental Ic{d) ob- 
served in (001) tilt grain boundary junctions can only be 
understood in terms of interface defects, regardless of the 
0P42 Of course, in Bi2212, it has long been established 
that a significant amount of static defects exist not only 
on the top cleaved surface;^ but throughout even the 
best single crystals i^^i^^*^ Such defects combined with 
significant oxygen loss at the grain boundaries may be 
the most important source of the observed 7r-junctions 
in the tricrystal and tetracrystal experiments, account- 
ing for the apparent inconsistency between those and the 
c-axis bicrystal and cross-whisker experiments4iiSii2iSi 

Finally, we recall that 7r-junctions are now rou- 
tinely produced using junctions of the SNS and 
SINIS varieties constructed from conventional 
superconductorSf^ii^Si^^i^^iSiS so that a close reex- 
amination of the tri- and tetracrystal experiments 



is warranted. This is especially true in the case of 
the hole-doped cuprates. We note that not only are 
TT-junctions readily obtained with Nb, Al and their 
oxidesj^ but the associated zero-bias conductance 
peaks and higher harmonic order Josephson phase 
relations in such conventional junctions have also been 
observed,a^2ia^lSia as in YBCO.^ 



V. SUMMARY 

We calculated the critical current Ic versus (001) in- 
plane tilt angle 9 for symmetric and asymmetric grain 
boundary junctions, assuming the order parameter has 
either the s-, extended-s-, or (i3,2_y2-wave form. We used 
two Fermi surfaces appropriate for hole-doped cuprates, 
one appropriate for electron-doped cuprates, and an ex- 
treme Fermi surface for comparison purposes, and took 
account of the surface boundary conditions appropriate 
for the interfaces. We studied both tunneling limits of 
specular and random tunneling. An important feature 
of our results is that the presence of the surface bound- 
ary conditions allows for a finite d-wave critical current 
across junctions in which the microscopic tunneling is 
random. This is especially important because the tun- 
neling across nearly all in-plane (001) tilt junctions is 
expected to be random on an atomic scale, even within 
each reconstructed facet. 

Although our results for electron-doped cuprates are in 
complete agreement with the interpretation of the obser- 
vation of the tricrystal experiments on those compounds, 
our results for hole-doped cuprates differ qualitatively 
with those obtained from Ginzburg-Landau models^iS 
and indicate that the interpretation of the tricrystal 
experiments is complicated, depending upon strong as- 
sumptions about the details of the Fermi surface topol- 
ogy and the nature of the tunneling processes. In ad- 
dition, corresponding problems in the interpretation of 
the tetracrystal experiments on hole-doped cuprates per- 
sist, leading to expectations of an lc{0) inconsistent with 
many experiments^ We urge further in-plane (001) tilt 
grain boundary experiments using films deposited by liq- 
uid phase epitaxy to be made, and for high resolution 
transmission electron microscopy studies of the pertinent 
grain boundaries to be presentedi^*^!^ 



Acknowledgments 

The authors would like to thank T. Claeson, S. E. 
Shafranjuk, and A. Yurgens for useful discussions. 



* Electronic address: |garnold@nd.edu] 



Electronic address: |kiemm@phys.ksu.edu| 



12 



^ C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000). 

^ H. Hilgenkamp and J. Mannhart, Rev. Mod. Phys. 74, 485 

(2002) . 

^ K. A. Miiller, Phil. Mag. Lett. 82, 279 (2002). 

* R. A. Klemm, Phil. Mag. 85, 801 (2005); Int. J. Mod. 

Phys. B 12, 2920 (1998). 
^ R. A. Klemm in Nonequilibrium Physics at Short Times 

Scales. Formation of Correlations, edited by K. Morawetz 

(Springer, Berlin, 2004), pp. 381-400. 
® P. Chaudhari and S. Y. Lin, Phys. Rev. Lett. 72, 1084 

(1994) . 

C. C. Tsuei, J. R. Kirtley, C. C. Chi, L. S. Yu-Jahnes, A. 

Gupta, T. Shaw, J. Z. Sun, and M. B. Ketchen, Phys. Rev. 

Lett. 73, 593 (1994). 
^ J. R. Kirtley, C. C. Tsuei, H. Raffy, Z. Z. Li, A. Gupta, J. 

Z. Sun, and S. Megtert, Europhys. Lett. 36, 707 (1996). 
^ C. C. Tsuei, J. R. Kirtley, G. Hammer, J. Mannhart, H. 

Raffy, and Z. Z. Li, Phys. Rev. Lett. 93, 187004 (2004). 
^° C. C. Tsuei and J. R. Kirtley, Phys. Rev. Lett. 85, 182 

(2000). 

" R. R. Schulz, B. Chesca, B. Goctz, C. W. Schneider, A. 
Schmall, H. Bielefeldt, H. Hilgenkamp, J. Mannhart, and 
C. C. Tsuei, Appl. Phys. Lett. 76, 912 (2000). 
B. Chesca, K. Ehrhardt, M. Mofile, R. Straub, D. Koelle, 
R. Kleiner, and A. Tsukada, Phys. Rev. Lett. 90, 057004 

(2003) . 

A. G. Sun, D. A. Gajcwski, M. B. Maple, and R. C. Dynes, 
Phys. Rev. Lett. 72, 2267 (1994); A. G. Sun, A. Truscott, 
A. S. Katz, R. G. Dynes, B. W. Veal, and G. Gu, Phys. 
Rev. B 54, 6734 (1996). 

R. Kleiner, A. S. Katz, A. G. Sun, R. Summer, D. A. 
Gajcwski, S. H. Han, S. I. Woods, E. Dantsker, B. Chen, 
K. Char, M. B. Maple, R. G. Dynes, and J. Clarke, Phys. 
Rev. Lett. 76, 2161 (1996). 

A. S. Katz, A. G. Sun, R. C. Dynes, and K. Char, Appl. 
Phys. Lett. 66, 105 (1995). 

1^ M. Mofile and R. Kleiner, Phys. Rev. B 59, 4486 (1999). 
S. I. Woods, A. S. Katz, T. L. Kirk, M. C. de Andradc, M. 

B. Maple, and R. C. Dynes, IEEE Trans. Appl. Supercond. 
9, 3917 (1999). 

^® Qiang Li, Y. N. Tsay, M. Suenaga, R. A. Klemm, G. D. 
Gu, and N. Koshizuka, Phys. Rev. Lett. 83, 4160 (1999). 
Y. Takano, T. Hatano, S. Kawakami, M. Ohmi, S. Ikcda, 
M. Nagao, K. Inomata, K. S. Yun, A. Ishii, A. Tariaka, 
T. Yamashita, and M. Tachiki, Physica C 408-410, 296 

(2004) . 

^° R. A. Klemm, Phys. Rev. B 67, 174509 (2003). 

Yu. I. Latyshev, A. P. Orlov, A. M. Nikitina, P. Monceau, 
and R. A. Klemm, Phys. Rev. B 70, 094517 (2004). 
M. Sigrist and T. M. Rice, J. Phys. Soc. Jpn. 61, 4283 
(1992). 

Y. Tanaka and S. Kashiwaya, Phys. Rev. Lett. 74, 3451 

(1995) ; Phys. Rev. B 53, R11957 (1996); ibid. 56, 892 
(1997). 

2" Y. Asano, Phys. Rev. B 64, 014511 (2001); Y. Asano, 
Y. Tanaka, and S. Kashiwaya, Phys. Rev. B 69, 134501 
(2004). 

Yu. S. Barasli, H. Burkhardt, and D. Rainer, Phys. Rev. 
Lett. 77, 4070 (1996). 
^® M. S. Kalenkov, M. Fogelstrom, and Yu. S. Barash, Phys. 
Rev. B 70, 184505 (2004). 

E. Il'ichcv, M. Grajcar, R. Hlubina, R. P. J. IJsselsteijn, H. 
E. Hoenig, H.-G. Meyer, A. Golubov, M. H. S. Amin, A. M. 



Zagoskin, A. N. Omclyanchouk, and M. Yu. Kupriyanov, 
Phys. Rev. Lett. 86, 5369 (2001). 

S. Shirai, H. Tsuchiura, Y. Asano, Y. Tanaka, J. Inouc, 
Y. Tanuma, and S. Kashiwaya, J. Phys. Soc. Jpn. 72, 
2299 (2003). In their Eq. (15), these authors defined their 
quasiparticle dispersion to he Cq = —{t + Ci)^ — Calg — 
(/It + 8Wn), with r/g = 2(cos(/a; + cosg^), their Eq. (17), 
in units of the lattice constant a = 1. Following their 
Eq. (21), the definitions t = W, (i = -O.Wt, (2 = O.Ot, 
II = — 0.2833t, and the hole density n = (X^o- = 0.85 
were given, where the summation is only over the two spin 
states. Putting these numbers into the above dispersion, 
one obtains tq — — 1.606t(cos + cosqy + 4.06), which 
never vanishes, and hence does not exhibit a Fermi sur- 
face. However, the q-independent part of this expression 
is inconsistent with their Eq. (4) for the effective chemical 
potential, jl = fi + W ^ {np,a}, where the summation is 
over the four near-neighbor in-plane sites denoted by p and 
the two spin states. Hence, the 8 in their Eq. (15) should be 
replaced with the number 4 to avoid overcounting the two 
spin states. With this modification, the correct dispersion 
of Shirai et al. becomes eq = — 1.606t(cos g^i-l-cos gj;-|-1.94), 
which has the Fermi surface at small pockets in the cor- 
ners of the first Brillouin zone shown in Fig. 2. Although 
Shirai et al. may have intended to include this factor of 4 
in their definition of n, which would have led to a Fermi 
surface resembling FS2 in Fig. 2, from the equations in 
their manuscript that they very likely used the above dis- 
persion corresponding in our notation to J|| = 201 meV, 
u = 0, and /tt = —1.94, where J|| can be determined from 
their value of A = 0. 07994, which we set equal to 10 meV. 
However, their strange results for the symmetric (mirror) 
and especially the straight (parallel) 18.4° junctions do not 
appear to correspond to any results by other authors for 
any FS studied. 

^® A. Bille, R. A. Klemm, and K. Scharnberg, Phys. Rev. B 
64, 174507 (2001). 

^° G. B. Arnold and R. A. Klemm, Phys. Rev. B 62, 661 
(2000). 

S. H. Liu and R. A. Klemm, Phys. Rev. Lett. 73, 1019 
(1994). 

32 X. Y. Cai, A. Gurevich, I.-F. Tsu, D. L. Kaiser, S. E. 
Babcock, and D. C. Larbalestier, Phys. Rev. B 57, 10951 

(1998). 

33 M. F. Chisholm and S. J. Pennycook, Nature (London) 
351, 47 (1991). 

^'^ X.-Z. Yan and C.-R. Hu, Phys. Rev. Lett. 83, 1656 (1999). 
3® G. B. Arnold and R. A. Klemm, unpublished. 
3® Y. Tanaka and S. Kashiwaya, Phys. Rev. B 58, R2948 
(1998). 

3'' R. A. Klemm, G. Arnold, C. T. Rieck, and K. Scharnberg, 

Phys. Rev. B 58, 14203 (1998). 
3* Yu. Eltsev, K. Nakao, Y. Yamada, I. Hirabayashi, Y. Ishi- 

maru, K. Tanabe, Y. Enomoto, J. Wen, and N. Koshizuka, 

Physica C 367, 24 (2002). 

V. Ambegaokar and A. Baratoff, Phys. Rev. Lett. 10, 486 
(1963); tbtd. 11, 104 (1963). 
*° Y. Ishimaru, J. Wen, N. Koshizuka, and Y. Enomot, Phys. 
Rev. B 55, 11851 (1997). 

X. F. Zhang, D. J. Miller and J. Talvacchio, J. Mater. Res. 
11, 2440 (1996). 

A. Gurevich and E. A. Pashitskii, Phys. Rev. B 57, 13878 
(1998) and references therein; ibid. 63, 139901(E) (2001). 



13 



K. M. Lang, V. Madhavan, J. E. Hoffman, E. W. Hudson, 
H. Eisaki, J. Z. Sun, and S. Mcgtcrt, Nature (London) 415, 
412 (2002). 

Y. Gao, P. Coppens, D. E. Cox, and A. R. Moodenbaugh, 
Acta Cryst. A 49, 141 (1993); P. Coppens, P. Lee, G. Yang, 
and H. S. Sheu, J. Phys. Chem. Solids 52, 1267 (1991); V. 
Petricek, Y. Gao, P. Lee, and P. Coppens, Phys. Rev. B 
42, 387 (1990). 

X. B. Kan and S. C. Moss, Acta Cryst. B 48, 122 (1992). 
D. Grebille, H. Leligny, A. Ruyter, P. Labbe, and D. 
Raveau, Acta Cryst. B 52, 628 (1996). 
H. J. H. Smilde, Ariando, D. H. A. Blank, G. J. Ger- 
ritsma, H. Hilgenkamp, and H. Rogalla, Phys. Rev. Lett. 
88, 057004 (2002). 

H. Hilgenkamp, Ariando, H.-J. H. Smildc, D. H. A. Blank, 
G. Rijnders, H. Rogalla, J. R. Kirtley, and C. C. Tsuei, 
Nature (London) 422, 50 (2003). 

A. Mathai, Y. Gim, R. C. Black, A. Amar, and F. C. 

Wellstood, Phys. Rev. Lett. 74, 4523 (1995); Anna Mathai, 

Ph. D. Thesis, University of Maryland, 1995. 

Y. Gim, A. Mathai, R. C. Black, A. Amar, and F. C. 

Wellstood, J. Phys. (Paris) I 6, 2299 (1996); Yonggyu Gim, 

Ph. D. Thesis, University of Maryland, 1996. 

J. J. A. Baselmans, A. F. Morpurgo, B. J. van Wees, and 



T. M. Klapwijk, Nature (London) 397, 43 (1999); J. ,J. 
A. Baselmans, B. J. van Wees, and T. M. Klapwijk, Appl. 
Phys. Lett. 79, 2940 (2001); J. J. A. Baselmans, T. T. 
Heikkila, B. J. van Wees, and T. M. Klapwijk, Phys. Rev. 
Lett. 89, 207002 (2002). 

R. Shaikhaidarov, A. F. Volkov, H. Takayanagi, V. T. Pe- 
trashov, and P. Delsing, Phys. Rev. B 62, R14649 (2000). 
J. Huang, F. Pierre, T. T. Heikkila, F. K. Wilhelm, and 
N. O. Birge, Phys. Rev. B 66, 020507(R) (2002). 
F. Giazotto, T. T. Heikkila, F. Taddei, R. Fazio, J. P. 
Pekola, and F. Beltram, Phys. Rev. Lett. 92, 137001 

(2004) . 

I. P. Nevirkovets, S. E. Shafranjuk, and J. B. Ketterson, 
Phys. Rev. B 68, 024514 (2003). 

T. M. Klapwijk, J. Supercond.: Incorp. Novel Magn. 17, 

593 (2004). 

T. Lindstrom, S. A. Charlebois, A. Ya. Tzalenchuk, Z. 
Ivanov, M. H. S. Amin, and A. M. Zagoskin, Phys. Rev. 
Lett. 90, 117002 (2003). 

T. Bauch, F. Lombardi, F. Tafuri, A. Barone, G. Rotoli, 
P. Delsing, and T. Claeson, Phys. Rev. Lett. 94, 087003 

(2005) . 



